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Abstract 

A model of an Ebola epidemic is developed with infected individuals structured 
according to disease age. The transmission of the infection is tracked by disease 
age through an initial incubation (exposed) phase, followed by an infectious phase 
with variable transmission infectiousness. The removal of infected individuals is 
dependent on disease age, with three types of removal rates: (1) removal due to hos¬ 
pitalization (isolation), (2) removal due to mortality separate from hospitalization, 
and (3) removal due to recovery separate from hospitalization. The model is applied 
to the Ebola epidemics in Sierra Leone and Guinea. Model simulations indicate that 
successive stages of increased and earlier hospitalization of cases have resulted in 
mitigation of the epidemics. 


1 Introduction 

The 2014-2015 Ebola epidemics in West Africa have been mathematically modeled in 
many studies with attention to various aspects of the disease epidemiology (References). 
In this study we focus attention on the variability of infectiousness of infected individuals 
throughout their disease course. It is recognized that the period of infectiousness follows 
a period of incubation lasting from 2 to 21 days, begins with the presentation of symp¬ 
toms, and that the level of infection transmission coincides with the severity of symptoms 
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Il34ll . Our analysis of infection transmission is based on the disease age of an infected in¬ 
dividual. Here disease age means the time since infection acquisition, which is assigned 
disease age 0. We assume that the disease age distribution of the level of infectiousness 
is a Gaussian distribution, rising gradually, peaking, and falling gradually for survivors. 
The value of tracking disease course by disease age is that removal of infected individu¬ 
als from the epidemic population can be tied to their level of infectiousness, and earlier 
removal can have significant benefit for mitigation of the epidemic. 

The Ebola epidemics in West Africa have developed with rapid increase in the cumu¬ 
lative number of cases from the spring of 2014 to the present. The most recent WHO data 
indicate that the epidemics may be subsiding li33l . In this study we analyze the dynamics 
of the epidemics in Sierra Leone and Guinea with a disease age model. Our goal is to 
connect the removal rates of infected individuals to the epidemic outcomes. Our simula¬ 
tions indicate that the mitigation of the epidemics is tied to an increased level and earlier 
hospitalization (isolation) of infectious individuals as the epidemics unfolded. This en¬ 
hanced hospitalization of identified cases can be tied to increased contact tracing, social 
awareness, and public health policies. The measure of these epidemics is tracked by the 
net reproductive values Ro derived from the models. Our simulations show that for both 
Sierra Leone and Guinea R 0 has decreased significantly to levels which indicate that the 
epidemics will subside in 2015. We fit our disease age structured model simulations to 
cumulative case data and extend our simulations forward in time to project the epidemic 
levels into 2015. 


2 The Ebola Model with Age of Infection 

The model consists of the epidemic population divided into susceptible S(t), infected 
/(f), and removed classes R(t ) at time f, with the key feature that infected individuals 
are structured by a continuous variable a corresponding to time since acquisition. The 
age of infection is divided into an early stage corresponding to an incubation period, and 
intermediate stage corresponding to an infectious period, and a late stage corresponding 
to a recovering period, as determined by the disease age of infected individuals. 

1. Initially there are So susceptible individuals with a small number of infected indi¬ 
viduals. The results are valid for any number of initial susceptibles with a given 
parameterization. 

2. The disease phases of an infected individual are governed by a Gaussian distri¬ 
bution of disease age a. The Gaussian distribution has mean a m , standard devia¬ 
tion a S d, and is multiplied by a scaling factor o 0 . It is assumed that the on-set of 
symptoms (coincident with infectiousness) is gradual after an incubation (exposed) 
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phase. The values a m and a sd should be viewed as effective average values for a 
typical infected individual. 

3. Infectious individuals are removed due to hospitalization/isolation at a rate Hi(a) 
per day after disease age a tl \- This rate takes into account time lags in identifying 
presentation of symptoms. 

4. Infectious individuals are removed due to unreported mortality (independent of hos¬ 
pitalization) at a rate p 2 (a) per day after disease age a^ 2 . This rate allows for delays 
resultant from improperly handled corpses. 

5. Infectious individuals are (effectively) removed due to recovery (independent of 
hospitalization) at a rate yU 3 (a) per day after disease age a M3 . The values of a m3 and 
a //2 are significantly greater than the value of a lt \. 

6. Infectious individuals are removed as a result of increased and earlier hospitaliza¬ 
tion at a rate yU 4 (a, t) per day dependent on disease age a and time t > t\ during 
the course of the epidemic. The time t\ is chosen based on case data indicating a 
lessening of the cumulative number of cases. We allow for a succession of times 
t\ < t 2 < ..., when additional increased and earlier hospitalization occurs during 
the course of the epidemic. 


The formulas for /ii(a), /i 2 (a), /r 3 (a), /x 4 (a,f) are 

fii(a) = 0 if a < ani, pi(a) = fi\o if a > a^i 

/i 2 (a) = 0 if a < a^ 2 , /r 2 (a) = /i 2 o if a > a^ 2 

^ 3 ( 0 ) = 0 if a < a M3 , /r 3 (a) = ^i 3 o if a > a^ 3 

/U 4 (a, t) — 0 if t < ti or a < a^ 4 , /U 4 (a, t) = /x 40 if t>t\ and a > a^ 4 . 
The probabilities of removal by disease age a for Hj, j = 1, 2, 3, 4, are 

PHj(a) = / Pj(b) exp ( — [ (ni(z) + fj, 2 (z) + p 3 (z))dz\db, j = 1, 2, 3, 


( 1 ) 

( 2 ) 

(3) 

(4) 


p/r 4 (a,f)= / /x 4 (&, f) exp 


(Hi(z) + n 2 (z) + // 3 (z) + n±{z,t))dz ) db, t > ti. 
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3 The Equations of the Model 


The equations of the model are based on the disease age density i(a, t ), where a is disease 
age and t is time. Infected individuals have disease age 0 at time of acquisition. The 
number of exposed (incubating) infected individuals at time t is E(t) = ["’ i(a,t)da, 
where a* is the approximate disease age at which the infectious phase begins. The number 
of infectious individuals at time t is I(t) = [('' i(a, t)da, where a r is the approximate 
disease age at which infectiousness ends. The number of susceptible individuals at time t 
is S (t). The model equations are 




a(a)i(a,t)da I S(t),t > 0 


(5) 


d d 

—z(a, t) + t— i(a, t) = —(/ii (a) + /i 2 (a) + / 13 (a) + /Tt(a, t), t > 0, a > 0 (6) 
at oa 


i(0,t) = 


a(a)i(a,t)da ) S(t),t > 0 


(7) 


d 

—R(t)— / (/ii(a) + /i 2 (a) + // 3 (a) + /i 4 (a, t))i(a, t)da, t > 0 ( 8 ) 

dt J q 

5(0) = So, i(a, 0) = * 0 (a),a > 0, R( 0) = 0 (9) 

An analysis of the equations Q-Q and the formulas below are given in [[30]]. The 
ultimate course of the epidemic satisfies lim 4 ^. 0 O /(f) = 0 and lini/^oo S(t) > 0. The 
reproduction number R 0 of the model before time t\ is 

Ro = S 0 J a(a)ex-p ( —J (fj,i(b) +/j, 2 {b) + 113 (b))db\da. (10) 

and after time t\ is 

R 0 = S(t 1 ) f a(a) exp ( — f (ni(b) + ji 2 (b) + n 3 (b) + /i 4 (6, ti))db ) da. (11) 


The total cumulative number of cases at time t, both hospitalized (reported) and non- 
hospitalized (unreported) is 

cumT(t) = Sq — S(t) = f i(0,s)ds. 

Jo 

The epidemic attack ratio at time t, 0 < t < t\, is 

AR(t) — [ i(0,t)dt/So- 
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The cumulative number of hospitalized cases at time t, 0 < t < fi, is 



The cumulative number of unreported cases at time t is cumT{t ) — cumH{t). The cumu¬ 
lative number of mortality cases (independent of hospitalization) at time t, 0 < t < ti, 
is 



The cumulative number of recovered cases (independent of hospitalization) at time t, o< 

t < ti, is 



Similar formulas can be obtained for AR(t ), cumH{t), cumM{t), and cumR{t ) for t > 
ti. 

4 Simulations of the Model for the Ebola Epidemic in 
Sierra Leone 

The model Q-Q is simulated for the 2014-2015 Ebola epidemic in Sierra Leone. For 
this simulation the epidemic is divided into 3 phases: May 27, 2014 to September 1, 
2014, September 1 to December 14, 2014, and December 14, 2014 forward in time. The 
simulation phase time values are t 1 = September 1, 2014 and t 2 = December 14, 2014, 
where each phase represents an increased removal of cases due to a greater number of 
cases hospitalized, with earlier hospitalization. The removal rates //1 (a), /j 2 (a), /r 3 (a) 
hold for the first, second, and third phases (t > 0). The removal rate // 4 (a, t) holds for the 
second and third phases ( t > /, ). Additionally, there is the removal rate firja. t ) holding 
for the third phase t >t 2 having the form 

/x 5 (a, t) = 0 if t < t 2 or a < a m5 , /i 5 (a, t ) = /i, 50 if t > t 2 and a > a m5 . (12) 

The parameters for the simulation are given in Table 1. The model parameters are 
chosen by fitting the model solution to cumulative reported case data from WHO situation 
reports Il33l . The graph of the disease age dependent transmission function a(a ) is given 
in Figure |T] The simulation of the model is illustrated in Figure [2j The phase times H and 
t 2 are chosen based on fitting this data to a solution of the logistic equation. The concavity 
change in the data corresponds to the transition over time of Rq from > 1 to < 1, which 
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corresponds in turn to the concavity change in the fitted logistic equation solution. The 
logistic equation and its solution are 

P'(t ) = r pm-P(t)/C), = (13) 

The dashed graph in Figure [2] is the solution of ([D) with r = 0.031, C = 11,650, and 
P(0) = 63. The dashed vertical line in Figured is the concavity change in P(t) at 
i m 169, Pit ) = C/2. In Figure [2] the simulation is projected forward in time to June 1, 
2015. The dependence of Rq during the first phase between May 27, 2014 and September 
1 , 2014 as a function of a M i (the earliest disease age for hospitalization in the first phase) 
and fi 10 (the hospitalization rate after disease age a /(l in the first phase) is illustrated in 
Figure [3] (all other parameters are as in Table 1). Figure [3] demonstrates the importance 
of smaller values of a ll: \ and larger values of fi \o for mitigation of the epidemic. The 
simulation for Sierra Leone is summarized in Table 2. 


5 Simulations of the Model for the Ebola Epidemic in 
Guinea 

The model (|5])-(|9]) is simulated for the 2014-2015 Ebola epidemic in Guinea. The param¬ 
eters for the model simulation are chosen as for Sierra Leone, based on WHO cumulative 
reported case data ll33ll . For this simulation the epidemic is divided into 3 phases: March 
26, 2014 to September 28, 2014, September 28, 2014 to December 21, 2014, and De¬ 
cember 21, 2014 forward in time. The simulation phase time values are t \ = September 
28, 2014 and to = December 21, 2014. The parameterization for the simulation is the 
following: S 0 = 10.6 x 10 6 ; to (a) is the same as for Sierra Leone (Table 1); the disease 
age dependent transmission rate a{a ) is the same as for Sierra Leone (Table 1) except 
a 0 = 0.49 x 10 -6 ; a M i, 2 , a m3 , a /l4 , a M5 are the same as for Sierra Leone (Table 1); the 
removal rate parameters are fii(a) = 0.35, /r 2 (a) = 0.22, /r 3 (a) = 0.30, /x 4 (a, t) = .053, 


yu 5 (a, t) = .045. The logistic equation parameters in (13) are r = 0.175, C = 3,752, 


and P( 0) = 65, with concavity change at t — 231, P(t) = C/2. The simulation of the 
model is illustrated in Figure |4j In Figure [4] the simulation is projected forward in time 
to June 1, 2015. The dependence of R 0 during the first phase between Mar 26, 2014 
and September 28, 2014 as function of i (the earliest disease age for hospitalization in 
the first phase) and /i 10 (the parameter for hospitalizationafter age disease a /( l in the first 
phase) is illustrated in Figure [5} The simulation for Guinea is summarized in Table 2. 
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6 Discussion 


We have presented a model of the 2014-2015 Ebola epidemics in Sierra Leone and Guinea 
structured by the disease age of infected individuals. The incorporation of disease age al¬ 
lows refinement of the disease transmission rate based on level of infectiousness. For 
Ebola the level of infectiousness is correlated to the level of symptoms, which begin 
gradually over a period of several days. The incorporation of disease age also allows 
refinement of removal rates due to hospitalization, mortality, and recovery correlated to 
disease age. Consequently, the significant advantage in hospitalization/isolation of in¬ 
fected individuals as soon as possible after presentation of symptoms is emphasized. 

Our model simulations are fitted to WHO cumulative reported case data ll33ll through 
February 22, 2015 and projected forward to June 1, 2015. The data indicate that the Ebola 
epidemics have subsided and will further subside to effective elimination in 2015 for these 
two countries. The reasons for the extinction of these Ebola epidemics are complex, but 
increased identification and isolation of infectious cases played a major role. Our model 
quantifies the significance of this increased identification/isolation. One way to quantify 
an epidemic is the epidemic net reproduction number R 0 , which is dynamic and evolves 
with changes in public health interventions and societal behavior. In our simulations 
we modeled this R 0 dynamic with two changes in hospitalization/isolation rates in two 
successive time periods. Our model simulations indicate that this removal of infectious 
individuals is key to elimination of the epidemic. This removal is quantified in our model 
in two ways: (1) with respect to the rate of infectious individuals hospitalized per day, and 

(2) with respect to an earlier disease age of the infectious individuals hospitalized. Both 
of these considerations can be influenced by public health policies and public awareness. 
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Figure 1: The disease age infection transmission function for Sierra Leone a (a) with 
a 0 = 1-27 x 10 -6 , a m = 16.5, a s d = 2.0, and the removal rate start ages a^i = 12, 
2 = 16, 3 = 17, 4 = 11, a^ 5 = 11. The infectious phase begins at ~ 10 days and 

ends at « 23 days. 
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Figure 2: Simulation of the Ebola epidemic in Sierra Leone with the parameters in 
Table 1. The dots are reported cumulative case data [f33l . The red graph is reported 
cumulative cases from the simulation. The blue graph is the total cumulative cases from 
the simulation (both reported and unreported). R 0 = 1.77 between May 27, 2014 and 
September 1, 2014. R 0 = 1.10 between September 1 and December 14, 2014. Rq = 0.72 
between December 14, 2014 and forward in time. 
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hospitalization/isolation parameter 



Figure 3: f? 0 for Sierra Leone during the first phase between May 27, 2014 and Septem¬ 
ber 1, 2014 as a function of a^i (the earliest disease age for hospitalization) and /ii 0 (the 
hospitalization parameter) with all other parameters as in Table 1. The blue surface is 
R 0 = , // io) and the red plane is R 0 = 1.0. The yellow dot corresponds to the 

Table 1 values a 7i] = 12 days, /i 10 = 0.35 per day, and R 0 = 1.77. 
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Figure 4: Simulation of the Ebola epidemic in Guinea. The dots are reported cumulative 
case data (331. The red graph is reported cumulative cases from the simulation. The blue 
graph is the total cumulative cases from the simulation (both reported and unreported). 
R 0 = 1.22 between March 26, 2014 and September 28, 2014. R 0 = 1.00 between 
September 28 and December 21, 2014. R 0 = 0.85 between December 21, 2014 and 
forward in time. 
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hospitalization/isolation parameter 



Figure 5: Ro for Guinea during the first phase between March 26, 2014 and September 
28, 2014 as a function of a M i (the earliest disease age for hospitalization) and // 1 0 (the 
hospitalization parameter) with all other parameters as in Figure [4} The blue surface is 
Ro = , // 1 o) and the red plane is R 0 = 1.0. The yellow dot corresponds to the 

values a M i = 12 days, /ii 0 = 0.35, and R 0 = 1.22. 
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Table 1: Parameters of the Model for the E 

:>ola Epidemic in Sierre Leone 

Parameter 

Parameter description 

Numerical value 

So 

initial susceptible population 

6 x 10 6 

i 0 (a) 

disease age density at time 0 

2.0 (0 < a < 9), 1.0 (9 < a < 23) 

h 

first improved reporting 

September 1, 2014 


second improved reporting 

December 14, 2014 

Q-ni 

disease age threshold of ^i(a) 

12 days 

M10 

/ti removal rate parameter 

0.35 per day 

Mi (a) 

hospitalization rate, £ > 0 

. . f 0.0, a < a„ 1 

v ' \ M10, a > a Ml 

a n2 

disease age threshold of /i 2 (0) 

16 days 

M20 

/t 2 removal rate parameter 

0.25 per day 

M2 (a) 

unreported mortality rate, £ > 0 

W (a) = | °'°’ “ ; a “ 2 

l M20, a > a M 2 

O/J 3 

disease age threshold of /x3(a) 

17 days 

M30 

/Z3 removal rate parameter 

0.35 per day 

Ms (a) 

unreported recovery rate, £ > 0 

Ms (a) = j °'°’ a J 0/13 
l M30, a > a^ 3 


disease age threshold of /r 4 (a, £) 

11 days 

M40 

/x 4 removal rate parameter 

0.13 

M4(^J 0 

increased hospitalization rate, £ > £ 1 

f 0.0, a < a^ 4 

M4(a, £) = < 0.0, t<t\ 

[ /i 4 Q, otherwise 

&/i5 

disease age threshold of /is (a, £) 

11 days 

M50 

/i 5 removal rate parameter 

0.13 


increased hospitalization rate, £ > £ 2 

f 0.0, (2 ^ 

/t 5 (a, £) = < 0.0, £ < £ 2 

[ /250, otherwise 

a (a) 

infection age transmission rate 

a(a) = a 0 Gaussian(a, a m , a s d) 

ao 

infection transmission scaling 

1.27 x 10~ 6 


infection transmission age mean 

16.5 days 

®sd 

infection transmission age sd 

2.0 days 
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Phase 

Sierra Leone 

Guinea 

First 

May 27, 2014 - September 1 

March 26, 2014 - September 28 

Ro 

1.77 

1.22 

total cases 

2,528 

1,742 

reported cases 

1,270 

1,234 

reported cases fraction 

0.50 

0.71 

Second 

September 1, 2014 - December 14 

September 28, 2014 - December 21 

Ro 

1.10 

1.00 

total cases 

14,976 

3,258 

reported cases 

8,108 

2,597 

reported cases fraction 

0.54 

0.80 

Third 

December 14, 2015 - February 22 

December 21, 2014 - February 22 

Ro 

0.72 

0.85 

total cases 

19,808 

3,573 

reported cases 

10,953 

2,864 

reported cases fraction 

0.55 

0.80 

Projection 

February 22, 2015 - June 1, 2015 

February 22, 2015 - June 1, 2015 

Ro 

0.72 

0.85 

total cases 

20,846 

3,983 

reported cases 

11,566 

3,211 

reported cases fraction 

0.55 

0.81 
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